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Mini-jet production in ultra-relativistic nuclear collisions leads to initial conditions 
characterized by large fluctuations of the local energy density (hot spots) and of 
the collective flow field (turbulence). Assuming that local equilibrium is reached on 
a small time scale, ~ 0.5 fm/c, the transverse evolution of those initial conditions 
is computed using hydrodynamics. We find that a new class of collective flow 
phenomena (hadronic volcanoes) could arise under such conditions. This could 
be observable via enhanced azimuthal fluctuations of the transverse energy flow, 
d 2 E ± /d<j>dy. 



1 Introduction 



At energies ji/i > 100 AGeV, mini-jet production in central nuclear collisions 
is expected El to be the primary dynamical source of a plasma of gluons with 
an initial energy density an order of magnitude above the deconfinement and 
chiral symmetry scale, e c ~ 1 GeV/fm 3 . Many observable consequences of the 
formation of this new phase of matter have been predicted based on a variety 
of assumptions □, and experiments are currently under construction to search 
for evidence for that so called quark-gluon plasma (QGP) at the Relativis- 
ts Heavy-Ion Collider (RHIC) at Brookhaven. Evidently, signatures depend 
sensitively on the assumed ensemble of initial conditions generated in such 
collisions. Most often it is assumed for simplicity that a homogeneous, cylin- 
drically symmetric, and longitudinally boost-invariant quark-gluon plasma is 
created and thus that signatures can be computed ignoring fluctuations of the 
initial conditions themselves. 

In this talk we point out, however, that the mini-jet formation mechanism 
leads to a rather inhomogeneous and turbulent ensemble of initial conditions 
that is characterized by a wide fluctuation spectrum of the local energy density 
(hot spots) and of the collective flow field (turbulence) . In this case, some of the 
proposed signatures will be washed out while new ones will certainly arise. We 
show below that in fact a new type of collective "shock" phenomena may occur 
under these conditions that could be readily observed as unusual transverse 
energy fluctuations. 

"Speaker at Int. Conf. Nucl. Phys., Wilderness/ George, SA, March 10-16, 1996 




1 



This topic is well suited to this symposium honoring the 60th birthday of 
Professor Greiner because he was the first to propose in 1974 the "shocking" 
idea that nuclear shock waves may be produced in nuclear x:ollisions and serve 
as an ideal probe of the equation of state of dense matteru. Since the exper- 
imental discovery of nuclear shocks in the 1 AGeV range in 1984 and in the 
AGS energy range last year ^ much theoretical and experimental effort has 
been devoted to this topic over a large energy range: Ei a b — 100 AMeV to 
100 ATeV. The present work focuses on a new class of shock phenomena that 
could occur in quark-gluon plasmas near the top end of that energy scale. 

To illustrate the novel type of collective flow patterns we have in mind, 
the hydrodynamic □ evolution of initial conditions consisting of three static 
cylindrical "hot spots" are shown in Fig. 1. 
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Figure 1: Hadronic volcanoes erupt due to shock formation in regions where expanding shells 
of matter evolved from inhomogeneous "hot spots" initial conditions (left panels) intersect. 
The signature of such volcanoes is the strong anisotropic azimuthal angle dependence of the 
transverse energy flow dEx /d<j> (right panels). Contours of constant energy density, e(x, y, r), 
at the breakup time Tf are shown in the middle panels. The top row shows the evolution 
from a symmetric initial configuration of hot spots. The bottom rows shows the evolution 
from a more realistic asymmetric configuration. 



''see talks by J. Symons, J. Stachel, and P. Braun-Munzinger in these proceedings. 
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An ideal gas equation of state (p = e/3) is assumed. Shock waves are seen 
to be created where the three rapidly expanding shells of matter intersect at 
six points. The collective acceleration of matter in those shock zones can be 
clearly seen as strong peaks in the transverse energy flow. For the symmetric 
configuration (Fig. la), a "mercedes" pattern emerges (Fig. lb), leading to 
three primary and three secondary shocks in Fig. lc. For the asymmetric case 
in Figs. ld,e the azimuthal angle dependence of the shock peaks is spread out 
and modified. The search for analogous enhanced anisotropic transverse energy 
fluctuations (or hadronic volcanoes jj) as a function of rapidity and azimuth is 
the subject of the following discussion. 

2 HIJING Initial Conditions 

Mini-jet production is thought to be the dominant mechanism controlling the 
initial conditions in A + A because the pQCD inclusive cross section for jets 
with moderate p± > pq ~ 1 — 2 GeV rises to the value (Jjetijpa) > 10 mb at 
RHIC energies. Compelling evidence for pQCD mini-jet dynamics has been 
observed in pp and pp reactions at collider energies 

Due to the nuclear geometry, the total number of mini-jet gluons produced 
in central ^4 + ^4 collisions is expected to be ~ A 4 / 3 ~ 10 3 times larger than in 
pp collisions. This simple geometric enhancement causes the rapidity density of 
mini-jet gluons to reach dN g /dy ~ 300 — 600 in Au+Au depending on the mini- 
jet transverse momentum cut-off scale po — 1 — 2 GeV/c. The hadronization of 
this mini-jet gluon plasma via the string fragmentation mechanism in HIJING 
approximately doubles the final hadron transverse energy distribution due to 
the pedestal or "string" effect. In the pre-hadronization phase of the evolution 
we therefore include a "soft" component in terms of a background gas of soft 
gluons normalized such that together with the mini-jet component the final 
transverse energy distribution predicted by HIJING is reproduced. 

The average initial energy density of the mini-jet gluon plasma can be 
estimated using the Bjorken formula: (e(r)) w (po/tqitR 2 ) dN g /dy (tq/t) ~ 
40 (to/t) GeV/fm 3 . This applies until the thermalization time is reached and 
work due to expansion must be considered. 

To see that the initial conditions are in fact highly inhomogeneous and are 
dominated by a few "hot spots" , we must calculate the local energy density 
e(r, x_l) instead of averaging over the transverse coordinate as in the Bjorken 
formula. The local energy density and transverse momentum density at z = 

C T.D. Lee first used this term for anisotropy that could arise due to surface instabilities. 
In our case the "volcanoes" are dynamically generated from specific inhomogeneous initial 
conditions. 
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(the y cm = frame) can be computed from the HIJING list of produced gluons 

as 



e(V, x x ) 
m(r, xj_) 





(x_l - x_l q (t)) . (1) 



1 + (rp± a ) 



The sum is over the produced gluons with transverse and longitudinal momen- 
tum (pj_ a , Pza = P_La sinhy a ). The longitudinal and transverse coordinates 
of the production, (z a — 0,xj_ a ), are determined from the coordinates of the 
binary nucleon-nucleon collision from which the gluons originate. The above 
formula takes into account the free streaming of gluons not only along the z 
direction via the volume element rAy but also in the transverse direction via 
x±a(t) = x ±Q + vj_ Q r, where v ±Q = p± a /p a . The factor 1/(1 + (h/Tp ±a ) 2 ) is 
the formation probabilityQ of the gluon in the comoving frame. High-pr gluons 
are produced first and gluons with lower pt later according to the uncertainty 
principle leading to the so-called inside-outside picture of the dynamics. Before 
r ~ h/p± a that component of the radiation field is still part of the coherent 
Weizsacker- Williams field of the passing nucleitO. 

The expression (|l|), when averaged over transverse coordinates, recovers 
the Bjorken expression (e(r)) sa £o( r o/ r ) with vanishing average transverse 
momentum density, (to(t)) = 0. 

To study the transverse coordinate dependence of the initial conditions, 
we must specify a transverse, Arj_, and longitudinal, Ay, resolution scale. 
The densities coarse-grained on that resolution scale arc obtained from ([!]) 
substituting 



6 2 (x ± - X±a (T))S(y a )^ V y ' ]yaU Q(Ar ± /2-\ X±a (T)-x ± \) . (2) 

Ar, Ay 



To determine the relevant resolution scale, we note that at time r the local 
horizon for any gluon in the comoving (y g = 0) frame has a radius cr. Thus 
at the thermalization time, r t h, each gluon can be influenced by only a small 
neighborhood of radius CT t h ~ 0.5 fm of the mini-jet plasma. We take the 
transverse resolution scale to be the maximal causally connected diameter, 
Ar± = 2cr t /j w 1 fm and the rapidity width to be Ay = 1. Gluons moving 
with larger relative rapidity y are produced later, T(y,pj_) ~ cosh(y)/px due 
to time dilation at the boundary of the local horizon. The above choice of the 
resolution scale is the most optimistic from the point of view of minimizing 
fluctuations of the thermodynamic variables between the causally disconnected 
domains in the transverse plane. We note that the number, Nd(r) = (R/t) 2 , 
of such causally disconnected domains in a nuclear area, irR 2 , is initially very 



4 



large. Even at later times, t <~ 3 fm/c, when the mean energy density falls 
below e c , several disconnected domains remain. 

The soft component is modeled by s=s 1700 low-pj^ gluons per unit ra- 
pidity with a Gaussian transverse momentum distribution with rms p± = 0.3 
GeV/c. This leads to about 500 GeV of transverse energy per unit rapidity 
as needed to reproduce the HIJING final dEr/dy. However, these soft gluons 
only add a relatively small, approximately homogeneous contribution to the 
energy density, e 8 w 4GeV/fm 3 , at r t h — 0.5 fm/c. 

3 Hot Spots and Turbulence 

Hot Spots and Turbulent Minijet Initial Conditions t=0.5 fm/c 




20 40 0.2 0.4 0.6 
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Figure 2: (a) Energy density distribution of hot spots in a typical HIJING event. Note the 
large fluctuations of the initial energy density across the transverse plane, (b) Momentum 
density fluctuations indicate considerable initial state turbulence, (c) Distribution of energy 
density on a Ar ± = 1 fm transverse coordinate resolution scale averaged over 200 HIJING 
events. Note that this resolution scale corresponds to the maximal allowed one by causality 
at proper time, t = 0.5 fm/c. (d) The thermal profile distribution on the same resolution 
scale has a rather large width. The mean and shape of course depend also on the equation 
of state (assumed here to be an ideal gas). 
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The inhomogeneous nature of the mini-jet plasma initial conditions in 
central Au + Au collisions at RHIC is illustrated in Fig. 2. In Fig. 2a the 
energy density profile at an assumed thermalization time r = 0.5 fm/c of a 
typical HIJING (b = 0) event at y/s — 200 AGeV is shown. The plasma in 
the central rapidity slice (Ay cm = 1) exhibits large fluctuations (hot spots) 
because the average number of hard mini-jets produced per nucleon is only 
~ 1 per unit rapidity at RHIC energies. The hot spots are also associated with 
a chaotic velocity field (turbulence) as seen in Fig. 2b. When averaged over 
many events the distribution of the proper energy density and temperature 
on a transverse resolution scale Ar± = 1 fm are shown in Figs. 2c, d. Note 
that at that formation time the plasma cannot be characterized by a unique 
temperature. In fact the widths of those distributions which are controlled 
by the Glauber geometry of finite nuclei and the size of the mini-jet cross 
section obviously cannot be neglected. Signatures of plasma formation will 
differ considerably in this turbulent gluon scenario than in the conventional 
homogeneous hot-glue scenario. For example high-p^ direct photon production 
and heavy-quark production will be greatly enhanced in the hot spots u. In 
the next section we find that hadronic volcanoes similar to those in Fig. 1 also 
arise. 

4 Hadronic Volcanoes and Transverse Energy Fluctuations 

In Figure 3 the azimuthal angle dependence of the transverse energy flow, 
dE±/dyd(j>, is shown for a typical HIJING initial condition. The left panels 
(Figs. 3a, c) correspond to an ideal gas equation of state, while the right panels 
(Figs. 3b, d) correspond to a first order Bag model equation of statea. The top 
two panels evolve the initial distribution in Fig. 2a with the momentum field 
in Fig. 2b artifically set to zero (m(xj_,ro) = 0). This static inhomogeneous 
initial condition is most similar to the one in Fig. 1. The initial isotropic nature 
of the thermally folded transverse energy distribution is seen by the dashed 
lines in Figs. 3a, b. The lower panels (Figs. 3c, d) correspond to the evolution 
of the turbulent initial condition in Figs. 2a, b. The turbulence of the initial 
condition is revealed by the anisotropy of the dashed curves. 

In all cases the magnitude of the transverse energy decreases with time 
due to the work done by the plasma undergoing longitudinal (Bjorken) expan- 
sion. More work is done by an ideal gas case than ajplasma with first order 
transition because the latter equation of state is softerti In the static hot spot 
case the production of hadronic volcanoes is indeed similar to Fig. 1. In the 
turbulent case, the effect is obscured somewhat by the large initial anisotropy. 
Nevertheless, in the turbulent case also several sharper peaks appear. Evidence 
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for the collective origin of the hadronic volcanoes in the turbulent case must 
be looked for in higher order Et correlation analysis, e.g., DSt{4>)Et{^)) ■ For 
a more detailed discussion of these results we refer to Ref.B. 
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Figure 3: The evolution of transverse energy cLEt / ! dyd4> at y = is shown for different 
inhomogcneous Hf JING initial conditions and equations of state. In parts (a,b) the fluid 
velocity field shown in Fig. 2b is set to zero while in parts (c,d) the turbulent velocity field is 
included. The evolution in parts (a,c) assumes an ideal gas law p = e/3, while in parts (b,d) 
a first order Bag model transition is included. The dashed curve is the initial transverse 
energy distribution at the onset of hydrodynamic evolution. In the static hot spot cases the 
initial transverse energy is uniformly distributed in azimuth, while in the turbulent case large 
initial fluctuations are generated by the mini-jets. The solid curves show the final transverse 

energy distribution. 
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